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Abstract 

We present a class of photonic lattices with an underlying symmetry given by a finite-dimensional 
representation of the 2-hlD Lorentz group. In order to construct such a finite-dimensional repre¬ 
sentation of a non-compact group, we have to design a T^T-symmetric optical structure. Thus, the 
array of coupled waveguides may keep or break 'PT-symmetry, leading to a device that behaves like 
an oscillator or directional amplifier, respectively. We show that the so-called linear PT-symmetric 
dimer belongs to this class of photonic lattices. 
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I. INTRODUCTION 


Photonic lattices provide a solid platform for the simulation of relativistic and quantum 
physics M. The use of finite or pseudo-infinite arrays, with an adequate set of effective 
refractive indices and coupling parameters, allows the optical realization of compact, SU(N), 
and non-compact groups, SU{1, 1), for example [5H7]. The ability to classify photonic lattices 
by their underlying symmetry opens the door to the simulation of a large class of quantum 
phenomena and the use of symmetries to design integrated photonic devices. 

Here, we are interested in the Lorentz group because it is ubiquitous in all linear theories 
where a state is represented as a wave function that belongs to a linear manifold with a scalar 
product, e.g., relativistic quantum mechanics, quantum mechanics, and electromagnetic 
theory in linear media [8HI3]. In optics, some matrix representations of the homogeneous 
Lorentz group have served as a bridge between the Jones and Mueller matrix calculus of 
polarization optics [IMH!, and have also been used in paraxial ray optics to describe multi¬ 
lens sytems, interferometers, laser cavities and multilayer systems [TH]. 

The unitary and irreducible representations of the Lorentz group are well known cni 
[la [13]. The Lorentz group is non-compact, thus, its unitary representations are infinite 
dimensional and we need infinite arrays of coupled lattices to create photonic lattices with 
this symmetry. In the laboratory, this means manufacturing an array with a large number 
of coupled waveguides and making sure propagated light never reaches one of the ends. The 
Lorentz group also has a family of finite-dimensional representations which are, however, 
non-unitary due to its non-compactness [121 EOj. The best way to construct these finite 
representations is to start from the su{2) Lie algebra with generators Jj and complexify it 
by defining Pj = iJj with j = 0,x,y. They all close a new Lie algebra that is isomorphic to 
the Lorentz algebra in 3-l-lD, so(3,1), containing various conjugated copies of the so(2,1) 
Lie algebra, {kQ,k^,ky} = {Jo,P^,Py}, {p,Py,Po}, {Jy,Po,P^}, that is, the Lorentz 
algebra in 2-l-lD. 

Taking into account that the propagation of electromagnetic fields through a photonic 
lattice with underlying SU{2) symmetry is modeled by the mode coupling matrix [7], 


HsU(2){z) — to Jo + A(J_ -|- j_|_) 
= coJq T 2ATj,, 


( 1 ) 

( 2 ) 


we find {k, Py, Pq} as the natural choice for an optical finite dimensional representation of 
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the 2+ID Lorentz group because it implies the change (u —?• ^7 in the effective refractive index 
parameter. This translates into waveguides with identical real part of the effective refractive 
index and non-identical imaginary part, that is, linear gain or loss at each waveguide [211E2]- 
In the following, we will study this optical finite representation of the Lorentz group. First, 
we are going to show that it belongs to the class of so-called T^T-symmetric optical systems. 
Then, we will provide a closed-form propagator for the device in terms of a response to 
impulse function. We will also show that the device behaves like an oscillator or directional 
amplifier for lattice parameters that keep or break the PT-symmetry, in that order, and, 
finally, provide a working example in the form of the so-called linear PT-symmetric dimer. 


II. PHOTONIC LATTICE SIMULATION 


Let us start with the coupled-mode differential set, 

-td,\Siz)) = H\S), (3) 


modeling the propagation of an electromagnetic field through an array of coupled waveguides. 
We have borrowed Dirac notation, such that column vectors are written as kets, 


\£{z)) = (£a(z),£i{z),...,£ 2 j(z)f , (4) 

= (5) 

/=0 

where the operation stands for transposition, and mode coupling matrices as operators 
illZ]. In particular, we are interested in waveguides with identical real part of the effective 
refractive index and linear gain or losses, that is, nonzero imaginary part of the effective 
refractive index, leading to the mode coupling matrix [H [ 2 T], | 22 ] , 


H 


Z7 {h- i)^-\ 


VJn{2j + 1 


-\-\Jh(2i -\-l- h)U'*' , 

+ A(J_ -|- J+), 


h) + 


( 6 ) 

(7) 


where the imaginary part of the effective refractive index and the effective coupling are given 
by the real parameters 7 and A, in that order, and we have used the SU{2) algebra that 
satisfies the commutation relations, [J+, j_] = 2 Jq and [Jq, J±\ = ± j±. Also, we have used 
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the step, h\j) = j|j), up-step, V^\j) = \j + 1), and down-step, V\j) = \j — 1) operators, that 
are diagonal, npper- and lower-diagonal square matrices of dimension 2 j, in that order. 

First, we want to emphasize that the effective refractive index of the n-th and {2j — n)-th 
waveguides are complex conjugates of each other, with j = 0,1/2,1,3/2 ,... and 2j -|- 1 is 
the number of waveguides in the device. This is a signature of the so-called PT-symmetric 
optical systems p3l - [28] . Thus, this photonic lattice has a dispersion relation, 

r 2 (m) = 0,0 {m- j), Oq = ^4A2 - 72 , ( 8 ) 


that may be real, 7 < 2A, complete degenerate, 7 = 2A, or purely imaginary, 7 > 2A. 

The impulse function of this array, that is, the field at the n-th waveguide given that the 
initial field impinged only at the m-th wavegnide, can be constrncted following the so-called 
Gilmore-Perelomov, symmetry based approach [7j, 


-Tn.n (^) 


= (n\e 




(9) 



fsA 


(^ 2 zAsin ^z^ 


[mj 

1 ) 



m+n 


X 


X iZnCos — 2;4-7sm — z] x 

V 2 ' 2 J 

( 4A2 .2^0 

X Km In; sm —z, 2j j , 


( 10 ) 


where the notation 



and Kn{x,p, N) stand for the binomial coefficient and Krawtchouk 


polynomials [22] • This impulse function is valid for any given parameter regime. 

A device mannfactnred with parameters that keep the PT-symmetry, 7 < 2 A, will show 
a 2t:/Oo periodicity, Fig. As we have said before, the fact that this is a non-nnitary 
representation translates into a photonic lattice with gain and losses. In practical terms, 
the PP-symmetric optical system behaves like a directional amplifier or attenuator. This is 
easily seen in Fig. ia) and Fig. [^d), where an initial field impinging the zeroth, m = 0 , or 
last wavegnide, m = 2j, will produce a peak in the intensity at the symmetric wavegnide, 
n = 2j — m, at a propagation distance Zp = {2k + l)7r/flo with A: = 0 , 1 , 2 ,..., 


^2j-m{zp) — ^2 ^2iAsin 2 ^ 


2 i 


m = 0 , 2j. 


( 11 ) 


In these particular cases it is simpler to see that for an input field impinging at the zeroth 
waveguide, m = 0 , the device behaves like an amplifier for half the propagation distance 
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FIG. 1. (Color Online) Light intensity propagation, for an initial field amplitude impinging 

just the (a) n = 0, (b)n = 3, (c) n = 5, (d) n = 10 waveguides of a lattice in the 'PT-symmetryc 
regime with parameters 7 = 0.4A and j = 5. 


and then like an attenuator for the rest of the propagation. Something similar occurs if the 
field impinges the last waveguide, m = 2 j, the device attenuates the input signal during 
the first half of the propagation length and then amplifies it during the second leg of the 
propagation. 
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FIG. 2 . (Color Online) Same as Fig for a completely degenerate lattice with parameters 7 = 2A 
and j = 5. 



In the fully degenerate case, 7 = 2A, it is straightforward to take the limit of the involved 
functions to find that the lattice looses its periodicity. 
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(12) 


X (ll; \^z^, 2j'j . 

Here, the device will work as a directional coupler, concentrating most of the field amplitude 
from the zeroth to j-th lattices if enough propagation length is given. Fig. In the broken 
T^T-symmetry region, 7 > 2 A, the directional coupling with amplification just becomes 
more pronounced. Fig. Note that in all figures the brightest value corresponds to the 
maximum intensity in the array, which is not limited to the unit value due to the presence of 
amplification. All figures were calculated via both a numeric solution of the mode coupling 
differential set and our analytic impulse function to good agreement. 
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FIG. 3. (Color Online) Same as Figfor a completely degenerate lattice with parameters 7 = 2.5A 
and j = 5. 

III. Pr-SYMMETRIC LINEAL DIMER 

Now, we want to work out a practical example and consider the case with j = 1 / 2 , 

- id^Soiz) = -^'ySoiz) + XSi{z), (13) 

-id^Si{z) = ^^£i{z) + XSoiz). (14) 

These equations describe the linear PT-symmetric dimer m [26], therefore, it possesses 
an underlying so( 2 , 1 ) symmetry realized non-unitarily Note that we can uncouple the 
dynamics via a second derivative with respect to the propagation distance, 

dl£j{z) + ^£tlSj{z) = 0, j = 0,1, (15) 
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and discover three different behaviors. In the PT-symmetric region, A > 7/2, each field 
amplitude behaves like an oscillator with frequency VLq/2, then, for lattice parameters A = 
7/2, they exhibit free particle-like behavior and, finally, the field amplitudes behave like 
a repulsive oscillator in the broken symmetry regime, A < 7/2. The response to initial 
impulses are given by Eq.(9) and can be simplified in the following form, 

( 16 ) 

( 17 ) 

( 18 ) 


To,o(^) = cos —z — sm —z, 
Z ^ Zq z 

Xo,i(^) =Xi^o{z), 

i2A 

^ / \ ^0 T • ^0 

Xi 1 z) = cos —z — — sm —z. 
’ ^ ’ 2 Qo 2 


(19) 


In order to visualize propagation in the linear PT-symmetric dimer, we will define a 
re-normalized light intensity at the n-th waveguide given an initial field impinging just the 
m-th waveguide in terms of the impulse function. 






( 20 ) 


The expression is the total light intensity at z and varies with the propagation 

distance due to the non-unitary propagation. Figure shows the re-normalized intensity, 
the left and right columns show the propagation of a light field impinging at the zeroth and 
first waveguides, in that order. The first row. Fig. ia) and|^b), shows a device in the VT- 
symmetric regime, 7 < 2A, where each field amplitude can be taken as the optical simulation 
of an harmonic oscillator with periodicity 27r/r2o; note that the re-normalized intensity is 
periodic but not harmonic. The second row. Fig. ic) and l^d), corresponds to the fully 
degenerate case, 7 = 2A, where the field amplitudes simulate free-particle propagation. In 
this case, the device asymptotically balances the re-normalized intensity at each waveguide 
for any given gain and coupling parameters. Finally, the third row. Fig. ie) and l^f), 
shows the case where the field amplitudes answer to an inverted oscillator dynamics where 
the scaled intensity have a more complex asymptotic behavior, 

( 21 ) 


lim VQfi{z) = lim 

Z^OO ’ z—)-00 ’ 

_ 7 + l^ol 
27 ’ 

Jim P„,i(^) = j™ Pm(2). 

^ 7 - IOqI 

27 


( 22 ) 

(23) 

(24) 
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These results help us recover the completely degenerate case, 


lim Vm,n{^) = 77 , 7 = 2 A. 

2;—)^00 5 ^ / 2 


(25) 



FIG. 4. (Color Online) Re-normalized intensity, Vm^niz), for a linear T^T-symmetric dimer with 
(a)-(b) 7 = 0.4A, (c)-(d) 7 = 2A and (e)-(f) 7 = 2.5A. The left and right columns show propagation 
for initial fields impinging at the m = 0 and m = 1 waveguides, in that order. Solid black lines show 
the re-normalized intensity at the n — 0 waveguide and blue dashed lines at the n — 1 waveguide. 


IV. CONCLUSION 

In summary, we have shown an optical finite-dimensional representation of the Lorentz 
group in 2 -l-lD. This is a non-compact group, thus, the representation is non-unitary and 
leads to an array of coupled waveguides with effective linear gain and loses, that is, the 
imaginary part of the effective refractive index is nonzero. The resulting photonic lattice is 
part of the so-called VT -symmetric optical systems and has two distinctive regimes where its 
dispersion relation is real and pure imaginary. Nevertheless, it is possible to use a symmetry 
based approach to construct a closed-form analytic impulse function that can be used to 
propagate any linear combinations of input fields in any given regime. We showed that in the 
regime that keeps the symmetry, with real dispersion relations, the device behaves like an 


























optical oscillator with amplification and attenuation. In the broken T^T-symmetry regime, 
the optical simulator delivers a directional coupler with amplification. 

As a practical example, we showed that the well known linear T^T-symmetric dimer 
has an underlying so(2,1) symmetry. Furthermore, we demonstrated that light propagating 
through it can act as an optical simulator of an harmonic oscillator, free-particle propagation 
or inverted oscillator in the T^T-symmetric, fully degenerate and broken symmetry regimes, 
in that order. Furthermore, the re-normalized intensity has a well behaved asymptotic 
behavior in the fully degenerate and broken symmetry regimes. 
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